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ABSTRACT 



Protostars are actively accreting matter and they drive spectacular, dynamic outflows, which evolve on timescales of years. X-ray 
emission from these jets has been detected only in a few cases and little is known about its time evolution. We present a new Chandra 
observation of LI 551 IRS 5's jet in the context of all available X-ray data of this object. Specifically, we perform a spatially resolved 
spectral analysis of the X-ray emission and find that (a) the total X-ray luminosity is constant over almost one decade, (b) the majority 
of the X-rays appear to be always located close to the driving source, (c) there is a clear trend in the photon energy as a function of the 
distance to the driving source indicating that the plasma is cooler at larger distances and (d) the X-ray emission is located in a small 
volume which is unresolved perpendicular to the jet axis by Chandra. 

A comparison of our X-ray data of the L1551 IRS 5 jet both with models as well as X-ray observations of other protostellar jets shows 
that a base/standing shock is a likely and plausible explanation for the apparent constancy of the observed X-ray emission. Internal 
shocks are also consistent with the observed morphology if the supply of jet material by the ejection of new blobs is sufficiently 
constant. We conclude that the study of the X-ray emission of protostellar jet sources allows us to diagnose the innermost regions 
close to the acceleration region of the outflows. 
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1- Introduction 



During the eariy stages of star formation the protostar is deeply 
embedded and therefore usually invisible at optical wavelengths; 
only infrared and radio emission, and potentially hard X-rays, 
penetrate the dense circumstellar environment. Yet, outflows es- 
cape the protostellar envelope and are detectable at various en- 
ergy bands, thereby announcing the birth of a new star At later 
stages of stellar evolution, when the accretion proceeds from a 
disk and the star becomes visible in the optical, outflow activ- 
ity is still observed. However, the driving mechanism of these 
outflows remains elusive; neither the acceleration nor the coUi- 
mation of the outflow are currently fully understood, magneto- 
centrifugally launched disk wind s, with a possible ste llar contri- 
bution are currently debated (e.g. lFerreira et al 

The most spectacular manifestations of these outflows are 
the condensations/shocked-regions termed Herbig-Haro (HH) 
objects. The proper motion of these knots within the outflows 
is generally in the range of a few 100km s '. Protostellar jets 
are intrinsically dynamic objects; their evolution is observable 
on timescales of a few years and models with variable ejection 
velo cities can su ccessfully explain some features of these jets 
(e.g. lRaga et al.ll201Ql . In these models, the overtaking of small 
slow blobs by faster more recently emitted blobs leads to shock 
fronts with shock velocities on the order of the amplitude of the 
velocity difference. 

HH objects with X-ray emission are a recently discovered 
phenomenon. About ten such objects have so far been de- 
tected among the hundreds of known HH objects; shock ve- 
locities around 500 kms"' are required to heat material to X- 
ray emitting temperatures in terms of simple shock models. 
HH 2 jPr avdo et al. 20 0l and die je t of L1551 IRS 5 (HH 154, 
iFavata et al.n2002 : .Ballv et alJl2003h marked the starting point 



of the X-ray discoveries. These outflows are driven by deeply 
embedded protostars (or their accretion disks). However, X-rays 
from the outflows are observed also from more evolved objects. 
The single classical T Tauri star DG Tau shows a complex X- 
ray morphology: There is the outer X-ray jet emission com- 
plex, the inner X-ray emission region locate d at a few 10 AU 
from the star and the stellar X- rav emission (j c iudel et al.ll2008l : 



ISchneider & Schmittll2008HGunther et ani2009ir 



In order to achieve high shock velocities within the jets, 
strongly varying outflow velocities are required to reach X -ray 
emitting temperatures. The models by iBonito et al.l (1201 Obh re- 
quire ejection velocities of a few 1000 km s ' in order to be rec- 
onciled with the X-ray observations. Such high velocities have 
not been detected in UV, optical or IR observations. As observ- 
able knots must experience at least one internal shock, these 
models predicted proper motion of the emission region with only 
a fraction of the initial flow velocity and are therefore also com- 
patible with the observations in the UV, optical and IR. Also, 
only a fraction of the total mass-loss carried by the outflow is 
probably required to explain the observed X-rays. Therefore, 
part of the high velocity material might have escaped detection. 
It is not clear, if the low number of detections is caused by the 
low X-ray surface brightness of these jets or if only a few out- 
flows actually emit X-rays at all. 

With our third epoch high-resolution Chandra observation 
of the L1551 IRS 5 jet (HH 154) we aim to determine the time- 
variability of the X-ray emission in order to constrain the origin 
of the X-ray emission. Our article is structured as follows: We 
introduce in the L155 1 IRS 5 region in the next section. In sect. [3] 
the observations and the data analysis are described. We then 
proceed to our results (sect. 2), which are discussed in sect. |5] 
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Fig. 1. Sketch of the region around the L1551 IRS 5 sources 
showing only a few components of the system (scales are only 
approximately preserved). The light gray regions indicate the 
Jets. The inner two [Fe ii] emission knots (PHKl, PHK2) from 
iPvo et al . (2009) are marked, the circumbinary disk is shown in 
dark gray. Note that the northern optical jet might actually be 
driven by the southern binary component (only one counter-jet 
is shown). The one arcsec wide stripes are used later in the article 
(sect.l4.2li and are shown here for reference. 



We review the implications of these results on current models in 
sect. |6] and close with a summary in sect. [7] 

As proposed by Bonit o etaP (l20I0bl) throughout the text 
the term "knot" describes a region of enhanced emission while 
"blob" refers to a moving gas clump which is not yet shocked, 
i.e. not observable in X-rays. 



2. L1551 IRS 5: Overview of previous observations 

The Lynds 1551 (L1551) star forming region (lLvndsllI962l) IS 
located at the southern end of the Taurus region at a distance 
of approximately 140 pc; within this region, a number of pro- 
tostellar objects and associated outflows have been found (e.g. 
iHavashi & Pvoll2009h . Fig. [1] shows a sketch of the immediate 
region around IRS 5. 

2.1. The L1 551 IRS 5 sources 

The term HH 154 describes the jet emanating from the sources 
collectively called L155 1 IRS 5 CStrom etaI..1976} . The infrai-ed 
source IRS 5 consists at least of a protostellar binary system and 
additional components are possible (see, e.g., the VLA data of 
[Rodriguez et al. 1998). Each core is surrounded by its own cir- 
cumstellar disk (~ 10 AU) and the co mplex is again embe dded 
within a large envelope (~ 10000 AU: iFridlund et al.ll2002h . All 
central sources are hidden by substantial absorption of at least 
Ay > 20 mag, but probably as much as 150 mag (e.g. ISnell et aP 
1198 5: Campbell et al. 1988). Therefore, the masses of the pro- 
tostars are uncertain, but fr om spectral rn apping of the reflected 
light escaping the envelope iLiseau et al.l (i2005i) derived masses 
of 0.3 and 0.8 Mq. These masses are consistent with the total 
system mass of 1.2 Mq derived bv Rodrig uez et all (120031) . The 
separation of the two sources, and consequently their jets, is only 
0.35" (50 AU). We will use the term IRS 5 for both sources 
jointly although most of the time the source of the northern jet is 
considered. 



2.2. The L1551 IRS 5 jet (HH 154) 

A double-lobed CO structure around LI 551 IRS 5 was first 
detected by Snell et al. (1980). Subsequent high resolution 
optical and near-infrared observations revealed two separate, 
westwards directed jets emanating from the immediate region 



around the two VL A sources. They can be t r aced out to abou t 
3" (420 AU, e.g. 'Fridlund & Liseau' '1998; 'itoh et al.' ' 20001) 
where the y become indistinguishable (Fridlund & Liseaul l998l: 
IPvo et al][2003h . The inner most part is on ly observable at ra- 



dio wavelengths (e.g. lRodrfguez et al]|2003l) . Since the northern 
jet is brighter and faster, it is believed to be responsible for the 
Herbig-Haro (HH) objects further downstream at distances be- 
tween a few and 1 2" (1700 AU). The j et inclination has been 
estimated to ~ 45° (IFridlund et al.ll2005b . 

At a distance of about 3 arcsec away from IRS 5 proper 
motion measurements of individual knots have been carried 
out on a baseline of 30 years, revealing substantial mo- 
tion of individual knots; the inner knots sho w the highest 
projected space-velocitie s of up to 300 km/s (IFridlund et al.l 
l200llBonito"erani200^ These values are approximately con- 
sistent with h ighest (projected) b lue shifted emission of up 
to 430 km/s (Fri dlund et all l2005b . However, high-resolution 
near-infrared [Fe II] 1.644 jum observations of the inner 
part showed an emissi on complex whi c h is v i rtually con- 
stant over four years ("P vo et al.l l2003l l2005l \2QQ% . The 
position-velocity diagrams (PVDs) show that the low-velocity 
component (LVC, v < 200 km s"') dominates the emission 
out to almost two arcsec, where the high-velocity compo- 



nent (HVC, 200 km s- 



< 450km/s) becomes domi- 



nant. High-resolution Hubble Space Telescope images taken 
in small band-pass filters (e.g.. Ha, [S ii]) can be explained 
by a light jet (i.e., less dense than t he ambient medium), hit- 
ting into a de nser ambient medium ( IFridlund &Liseaul 1 19981: 
Hart igan et al. i 2000). Spectroscopically, the outer knots show 
a line-width of 110 kms"', densities from a few 10"^ cm"^ to 
8 X lO-' cm'-^ and an excita tion rising with decreasing distance to 
IRS 5 (ILiseau et alJl200l and references therein). Concerning 
the nomenclature, the visual knots are designated F, E an d D (in 
increasing distance from IRS 5, see IFridlund et al1l200l. the in- 
ner near-infrared knots are termed PHK 1. . . 3 (Pyo et al. '200j 
in increasing distance to IRS 5 (cf. Fig. [TJ. Knot D coincides 
with PHK 3 and knot F with PHK 2. Whether IRS 5 is also driv- 
ing HH 28 and HH 29, which are located further downstream of 
HH 1 54, is not yet clear ( Devine et al. 1 999), possibly L 1 55 1 NE 
is their driving source. 

Throughout the text we will use the term HH 154 for all out- 
flow parts associated with the L1551 IRS 5 jet. 



2.3. X-rays from HH 154 

HH 15 4 was first discovered as an X-ray source by iFavata et al.l 
(I2OO2I) from an observation with XMM-Newton. Despite of the 
large PSF of XMM-Newton (~ 15 "), the authors coiTectly con- 
cluded from the spectral properties of the photons that the X-ray 
emission cannot be associated with IRS 5 and proposed knot D 
as the source of the X-ray emission; knot D is the brightest op- 
tical knot and probably the current terminal working surface. 
Higher resolution Chandra observations (see Fig. |2]l revealed 
that the X-ray source is actually located fu rther inward towa rds 
IRS 5 at a distance of 0.5-1.0" from IRS 5 (iBaUv et al.ll2003l) . A 
second epoch Chandra exposure showed a somewhat different 
morphology of the X-ray emission being more elongated than 
the 2001 ACIS data. This elongation has been interpreted as 
a moving X -ray knot with a p rojected space velocity of about 
330 km s ' (iFavata et al.ll2006l) . 

Based o n their initial detection of X-rays from HH 2 
iPravdo et al.l (I2OOII) proposed shocked, high-velocity knots as 
the explanation of the observed X-rays. The first analytica l de- 
scription of this process was presented bv lRaga et al.l(l2002l) and 
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Fig. 2. The ellipse for the global plasma properties, also, the location of the 2005 "knot" is shown. The crosses indicate the position 
of the radio sources and their size approximately indicates the position uncertainty. 



Table 1. Analysed X-ray observations of HH 154. 



Date 


Observatory 


Setup 


Obs-ID 


exp. time 


2000-09-09 


XMM-Newton 




0109060301 


56 ks 


2001-07-23 


Chandra 


ACIS-I 


1866 


80 ks 


2004-03-. . 


XMM-Newton 




0200810201... 


107 ks 


2005-10-2?" 


Chandra 


ACIS-I 


5381 


98 ks 


2009-12-29 


Chandra 


ACIS-S 


11016 


66 ks 



"This dataset consists of 1 1 short exposure distributed over six days 
in March 2004 where HH 154 is off-axis by about 2 arcmin. 



iBonito et alJ (l2004t) performed the first numerical hydrodynamic 
models with an emphasis on the X-ray emission. These models 
have been ex tended towards variable blob ejection velocities by 
iBonito et alJ (2010a b) and their analysis revealed that very fast 
blob velocities of more than 1000 km/s are needed in order to ex- 
plain the observations by shock heating, i.e., by internal shocks 
occurring when fast blobs overtake slower ones. An ejection "pe- 
riod" of two years matched the X-ray observations best. 



3. Observations and data analysis 

Table[T]lists all the available X-ray observations of HH 154. We 
used the ACIS-S detector for the third epoch Chandra exposure, 
since the back-illuminated ACIS-S chip has a higher sensitivity 
at lower energies than the front-illuminated ACIS-I CCDs. With 
the VFAINT mode, this setup provides a similar sensitivity as the 
longest ACIS-I exposure for the energies at hand; and an even 
higher sensitivity for plasma at cooler temperatures as expected 
for individual knots moving outwards and cooling. 

We used CIAO version 4.2 throughout the data analy- 
sis and followed the science threads on the CIAO webpagq^- 
The ACIS-S observation was reprocessed to account for the 
VFAINT-mod^ We experimented with pixel randomization, 
but since the relevant scales are usually at least twice the de- 
tector pixel scale, the effect of pixel randomization is virtually 
negligible. Therefore, we used the standard processing includ- 
ing pixel randomization. We explicitly note in the text where we 
expect this assumption to be invalid. 

In order to improve the astrometric accuracy, the three 
Chandra observations of HH 154 were aligned by calculating 
the centroids of the brightest sources detected by the CIAO tool 
celldetect and the photon events were reprojected so that the 
mean offset between these centroids vanishes. We use the off- 
set obtained from the three brightest sources weighted by the 



square-root of their count-number. These sources are also mem- 
bers of the Taurus star forming region and located on the ACIS- 
S part of the CCD array. Using more sources (detected with 
S/N>3) and equally weighted or photon number weighted means 
changes the offset by less than 0.5 pixel (0.25")- Thus, relative 
positions are at least accurate to within one pixel (0.49")- For 
a cross-check of the positions we calculated the centroid posi- 
tions of HH 154 using the photons within a circle of 2.5 pixel ra- 
dius centered on the brightest emission peak in the energy range 
0.5 - 3.0 keV. They coincide to within approximately 0.3", well 
within our estimated accuracy. Note that these centroids should 
coincide only in case of a stationary source. We give distances 
relative to the radio position with respect to the nominal po- 
sitions of the 2001 observation where the comparison sources 
show a good agreement with their optical positions but note that 
this position is accurate only within 0.5". 

We also retrieved the archival XMM-Newton data 
("obs^^idO 109060301) from 2000 for a spectral cross- 
check and eleven exposures during March 2004 
(obs-id 02008 10201... 02008 11 301), where HH 154 is lo- 
cated about two arcmin off-center for a luminosity check. 
SAS 9.(0 was used for the analysis of the XMM-Newton data. 
We extracted the source photons within a circle of 15" around 
the source position of HH 154 for the March 2004 exposures 
and derived the background from a nearby source-free region. 
We concentrated on the MOS data since the PN suffers high 
background levels (using standard filters only 40 ks PN on- 
time remain). However, both count-rates agree within their 
respective Icr ranges. We converted the count rate to luminosity 
by assuming the same spectral properties as during the 2001 
XMM-Newton observation. 



4. Results 

The X-ray images of HH 154 in the 0.5-3.0 keV energy band for 
the three available Chandra observations are shown in Fig.|2] To 
extract photons we use an ellipse with semi-axis lengths of 3.7" 
and 2.3", respectively, whose semi-major axis is aligned with 
the centroid in declination, i.e., which is aligned approximately 
with the jet axis and contains all photons attributable to the X- 
ray emission of HH 154. Using a nearby source free background 
regio n (no X-ray sources detected nor a 2MASS (Skrutskie et aT] 
120061) or SIMBAEQ source known), the expected background 
{E photon - 0.5 . . .3.0 keV) within the ellipse is 0.7 photons for 
the 2001 observation, 1.0 for the 2005 observation and 1.4 in the 
2009 observation. 



'http : // cxc . harvard . edu/ ciao/ 

^http : // cxc . harvard . edu/ ciao/threads/createL2/ 



'http : //xmiti . esa . int/ sas/ 

"^http : // simbad . u-strasbg . fr/ simbad/ 
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Fig. 3. Unabsorbed X-ray luminosity (0.5-lOkeV) of HH 154 
including the XMM-Newton data for the 2000 and 2004 data 
points. As the XMM-Newton data lacks sufficient spatial reso- 
lution to resolve the X-ray emission, the displayed luminosity 
pertains to the total observed X-ray emission for all datasets. 



4.1. Energetics 

We used XSPECvl2.5.0 (lArnaudlll99"6l) for the spectral mod- 
eling and assumed that the observed material can be described 
by op tically thin thermal plasma emission (APEC, ISmith et aP 
1200 lb and included absorption by neutral gas along the Une 
of sight in our fits. We set the abundance to half the solar 
value (using Anders & Grevesse 1989) since the XEST survey 
of the Taurus regio n found on average a sub-solar metallicity 
(iGudel et al.ll2007ah . 

Figure [3] shows the unabsorbed X-ray luminosity of HH 154 
during the last decade. It clearly indicates that the luminosity 
in the HH 154 region appears constant and deviations from the 
mean value do not exceed 22%. 

4.1 .1 . Global plasma properties 

The spectra extracted using the photons within the ellipse are 
shown in Fig. |4] however, we caution that the assumption of 
homogeneous plasma properties throughout the emitting region 
is probably not valid (see sect. 14.1.2b . Averaging over different 
plasma properties leads to an unstable fit with two solutions de- 
scribing the data reasonably well; the two possibilities are listed 
in Tab. |2] (top). One solution is only weakly absorbed and re- 
quires rather high plasma temperatures. Although this solution 
is statistically favored, we regard this solution as physically less 
plausible due to the following reasons. The Si lines at ~1.9 keV 
are not reproduced by this model but clearly present in the data 
(bottom panel of Fig. |4]i and the low Nh value contrasts the high 
absorption derived for the region close to the driving sources 
(Nh > 10^-^ ™d along the jet axis (A^^ > 8 x 10^' cm"^, 

e.g. Iltoh et alj|2000l: iFridlund et af]|2005h . In the following, we 
will therefore concentrate on the fit solution I with higher ab- 
sorption and lower plasma temperature. The plasma properties 
within the ellipse are compatible with each other for the individ- 
ual ex£osures(lo;}/rhe discrepancy between the values given 
bv iFavata etall (l2002l) and our values results partly from the 
change of metallicity (the respective 1 cr ranges overlap for solar 
metallicity). 

In order to estimate an upper limit for the presence of cooler 
plasma, we added a second temperature component to the fit and 




0.5 1 2 5 

Energy (keV) 




0.5 1 2 5 

Energy (keV) 



Fig. 4. Spectrum of all photons within the ellipse (see Fig. |2]l. 
Top: Individual spectra. Bottom: Co-added spectra (the Si-lines 
around 1.8 keV are clearly visible). 



fixed its temperature to 0.3 keV (0.2 keV). The Icr upper limit 
on the luminosity of this cool component is 1.4 x lO^^ergs ' 
(3.0 X 10^^ erg s"') when the absorption is forced to the value of 
the one temperature fit which also agrees with the optical/near- 
infrared value. Allowing the absorption to vary, the luminosity 
of the low temperature component decreases as the absorption 
also decreases for this two temperature component fit. 

4.1.2. Local plasma properties 

For a quantitative comparison of the three Chandra observations 
we divide the region around HH 154 into ~ 1" (2 ACIS pixels) 
wide spatial bins as indicated in Fig.|6] For a point source about 
80% of the photons are located in these 1" wide stripes. This 
procedure is essentially a projection of the photon number onto 
the flow-axes (i.e., the x-axis). Figure |6] shows the result for the 
individual exposures. Naturally, the exact values depend on the 
stripes used, therefore, we checked the results by shifting the 
stripes or using a different width of the stripes (about 90% of the 
photons of a point-source would be included in a three pixel wide 
stripe). Any property which depends crucially on the choice of 
the stripes is regarded as an unphysical artifact (Tab.|3)- 
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Table 2. Plasma properties of HH 154 with Icr errors. Spectra 
are binned to ~ 44 eV wide channels, i.e. about three times over- 
sampling the intrinsic energy resolution of the ACIS detectors, 
we used c-stat and the energy range 0.5-5.0 keV. Unabsorbed 
fluxes (0.5-10. keV) are given. 



Nh 


kT 


EM Lx 


(10^2 cm-2) 


(keV) 


(10^' cm-3) (10" erg/s) 





Composite spectrum (ellipse) 




Solution I 


1.1 ±0.1 0.6 ±0.1 


7.9±2.2 


82.7!^?:« 


obs2001^ 






89 7+^''-3 

"^•'-24.4 


obs2005>2l 






81-6lli 


obs2009>2l 




' -3.4 


77.6!^!:;; 


Solution II 


0.2 ±0.1 l.S+l\ 


1.9±0.2 


18 4+-' 




Individual spectra (ellipse) 




XMM (2000) 


1.0 ±0.1 0.5 ±0.1 


°-^-2.0 




obs2001 


1 1 6+° '' 

1.1 3 "■'-'-0.1 


8 9+'' * 


qn 9+37.1 
7;.z._33 2 


obs2005 


1.1 ±0.2 0.6±0.1 




79 8+2*'-'' 


obs2009 


1.0 ±0.2 o.6;!;-2 


7.0!^| 


73 4+32-7 




Composite spectra 






Eastern region 


l.4+°j 0.7±0.2 




78.1 


Western region 


O.e'ill 0.3±0.2 




5.0 



"Nh and fcr fixed to values from co-added spectra 



The mean energy of the photons in the individual strips along 
the flow axis is displayed in Fig. |5] (top panel), showing a clear 
decrease of photon energy with increasing distance to the driv- 
ing source(s). The slope of the mean energy depends on the de- 
tector's spectral response and on the stripes used. However, the 
spectral softening with increasing distance is independent of the 
detector response since only relative changes within the same 
detector are compared. 

For a cross check of the trend of the mean energy, we divided 
the emission region into an "eastern" and a "western" region 
so that the component close to the driving source(s) is associ- 
ated with the "eastern" region and the outer part of the emis- 
sion with the "western" region (these two regions essentially 
split stripes 2-7 into 2 separate regions); Table |2]lists the asso- 
ciated fit results. Since the total number of counts in the "west- 
ern" region is only 42 when summed over all ACIS exposures 
(^photon ~ 0.3-5.0 keV), we checked the results by fixing ei- 
ther the temperature or the absorption to the values obtained for 
flie left part of the emission, this results in a decrease of the Nh 
value consistent with Fig. |5] and the temperature for the right 
emission component, respectively. 

Figure |5] (bottom panel) shows the change required in either 
the plasma temperature or the absorption to explain the trend 
of the mean photon energy; fixing one parameter requires large 
changes of the other parameter (Nh = 2 x 10^^ — > cm"^ or 
/tr = 1.2 ^ 0.2 keV). 

The errors shown in Fig. |5](top) were obtained by simulating 
spectra containing a specified number of photons and then calcu- 
lating the mean energy range which contains 68 % of the trials. 
For photon numbers larger than ~20, the error (±0.1 keV) de- 
pends only very weakly on the number of photons; furthermore. 
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Fig. 5. Trend of mean energies, stripes are those of Fig. |6] 
Photons used within 0.5-3.0 keV. The red dash-dotted line is in- 
tended to guide the eye in the lower panel. Top: Thick sym- 
bols indicate data points with more than three photons while 
the smaller symbols refer to data points with fewer photons. 
Errors are obtained using the simulations of sect. 14.1.21 Bottom: 
Simulations for fixed temperature and fixed absorption along the 
jet axis. 



the error depends only weakly on the assumed plasma properties 
for the spectra at hand. 

4.2. Morphology along the jet axis 

The structure of the X-ray emitting region observed in 2009 nei- 
ther resembles the structure present in 2001 or 2005; it is more 
extended than the 2001 structure, but does not show excess emis- 
sion as far downstream as the 2005 exposure. The emission re- 
gion close to the driving source, which is also present in all 
previous exposures, is most notable. The new ACIS-S observa- 
tion does not show a clear knot westwards (downstream) of the 
main emission component as suggested by the 2005 image and 
there is no X-ray emission even further downstream as would 
be expected for a moving knot of constant luminosity. Note that 
the ACIS-S exposure is more sensitive to low energy photons 
than the 2005 ACIS-I exposure and that, according to Fig. |6] the 
photons soften with increasing distance to the driving source. 
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Table 3. Stripe properties (mean energy is in keV). 



Observation Stripe 







1 




2 




3 




4 


5 




6 




2001 


Counts 







17 




37 




4 


4 











Mean energy 






1.7 




1.4 




1.2 


1.1 




_ 




2005 


Counts 







22 




28 




9 


7 




4 






Mean energy 






1.6 




1.4 




1.2 


1.1 




1.2 




2009 


Counts 


3 




13 




26 




13 


8 




1 






Mean energy 1.2 




1.7 




1.4 




1.4 


0.9 




1.7 












Stripes shifted by 0.5" 






2001 


Counts 




7 




33 




16 




5 


1 




1 




Mean energy 




1.7 




1.5 




1.3 




1.2 


1.2 




1.0 


2005 


Counts 




4 




36 




14 




6 


7 




3 




Mean energy 




1.8 




1.5 




1.4 




1.2 


1.0 




1.6 


2009 


Counts 




3 




27 




19 




12 


1 




1 




Mean energy 




1.2 




1.6 




1.3 




1.1 


1.2 




1.7 



Table 4. Spectral models for sect. 14.2.11 



Model 


Nh 


kT 




(10^2 


(keV) 


Best tit values 


1.1 


0.6 


Fixed temperature 


0.0-2.0 


0.5 


Fixed absoiption 


0.6 


0.3-1.2 



Therefore, any emission with comparable properties as the pho- 
tons attributed to the "knot" should be detectable with the 2009 
ACIS-S observation. 

4.2.1 . Comparison of the Chandra observations 

A direct comparison of the photon numbers is not possible due to 
(a) different exposure times and (b) different detector responses. 
Nevertheless, we show in Fig. |6]the mean photon number ig- 
noring the differences in the sensitivity of the individual expo- 
sures. We note that all observations are compatible with this 
rou gh mean value using the error obtained by Gehrels weight- 
ing iG^dlll 9 8 6 ) . 

For a more detailed comparison of the observations, a spec- 
tral model is constructed (not fitted) in each stripe. We exper- 
imented with models which use the overall best fit values and 
with models which reproduce the trend in the mean energies (see 
Fig.|5]), they are listed in Tab.|4] As the predicted count numbers 
differ by less than one count, the statistical error overwhelms 
the error due to the unknown spectrum. Due to the different en- 
ergy response of the 2001 ACIS-I and 2009 ACIS-S, the scal- 
ing factors relative to the 2005 observation are 0.83. . . 0.87 and 
0.91. . . 1.20 depending on the assumed spectra (high photon en- 
ergy and low photon energy, respectively). 

The individual stiipe models are normalized so that the to- 
tal count number summed over the three ACIS observations in 
each stripe is conserved. Thus, the predicted total count number 
in a single stripe matches the observed value, which is statisti- 
cally the best estimate for the model normalization for constant 
emission. Figure |6] exemplarily shows the result for the model 
with kT = 0.5 keV and variable absorption along the flow-axis 
(which is virtually indistinguishable from the model with con- 
stant absorption and variable temperature). 

4.2.2. Time variable emission? 

To test if the observations are statistically consistent with the 
hypothesis of constant emission, we perform Monte-Carlo sim- 



ulations to estimate a confidence interval to accept or reject this 
hypothesis, i.e., to check if the observed photon distribution is 
an exceptional realization for a time constant emitting region. 
Since neither the time when new blobs appear nor their speed 
is predicted by cuiTent theories, the location of a new blob is 
not known a priori. Therefore, the statistical significance of any 
count number enhancement in a given region depends on the 
number of independent regions in which such an enhancement 
would be considered a knot. Essentially any stripe in Fig. |6]can 
be regarded as a possible region for a knot, however, the result 
does not depend strongly on, e.g., the inclusion of stiipes 1 or 
7, although the general result does depend strongly on the set of 
stripes used. To ensure that our results are not biased by the par- 
ticular selection of stripes, we repeated the simulation for differ- 
ent sets of stripes, which are two or three pixels wide and which 
have mutual offsets of one pixel. 

One simulation involves a set of three new exposures simu- 
lated as new individual photon numbers for each stripe. These 
simulated photon numbers were based on the expected photon 
number in that stiipe (sect. l4.2Tt and Poisson statistics. The in- 
dividual likelihoods, i.e., the likelihoods to observe exactly the 
simulated photon number in a given stripe-observation combi- 
nation, depend on the assumed expectation value for the pho- 
ton number in this particular stiipe. We derived this expectation 
value a posteriori from the simulated counts in a particular stiipe 
using the relative efficiencies from sect. 14.2.11 The total like- 
lihood for each simulation was then calculated as the product 
of the individual likelihoods. Thus, the fraction of realizations 
with a total likelihood better than the observed one can be in- 
terpreted as the probability for time variability. This approach 
does not include background events (which is minor effect on 
the ~ I % level) or other detector effects like alignment errors 
between the individual exposures. We find probabilities between 
approximately 50% and 96% depending on the exact stripes 
used. As the spatial distribution of the luminosity along the jet is 
unknown, we cannot decide if the set of stripes indicating time 
constant emission or the set of stripes indicating variable emis- 
sion, matches the real jet better Thus, we conclude that time- 
variable and time-constant emission are statistically acceptable. 

A hypothetical moving, fading and X-ray emitting knot 
would result in a low photon number at the western end of the 
box during the 2001 observations since at that time the knot was 
located more closely to the driving source, a higher photon num- 
ber during the 2005 observation, and virtually no photons by 
2009 due to cooling and expansion. The "knot" region (as indi- 
cated in Fig. |2| shows exactly this count number pattern. Thus, 
such knots are also compatible with the observations. 

In summary, the Chandra X-ray observations can be re- 
garded as statistically compatible with time constant X-ray emis- 
sion. However, the significance is poorly constrained as it de- 
pends crucially on the stripes used. The inability to find a clear 
sign for time dependence independent of the stripes used might 
be caused by the low number of counts since the 90 % confi- 
dence level easily covers a range almost twice as large as the 
value itself in the outer regions of the outflow. Therefore small 
scale time variability may be present, however, statistically time- 
dependent emission is not required. 

4.3. Extent perpendicular to the flow direction 

In order to check whether the Chandra exposures show evidence 
for an extent of the emission perpendicular to the jet axis, we 
have to take into account the exact position angle (PA) of the 
jet. Otherwise, the distribution of the photons around the jet axis 
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Fig. 6. One example of stripes ( 1 "width) used to extract count numbers, mean energies and spectra. In the bottom row, the projection 
onto the outflow direction (here assumed to coincide with the x-axis) is shown, the thick black line is the measured photon number, 
the thin blue dashed line indicates the mean photon number ignoring differences in the efficiency of the individual exposures. The 
red bars indicate the expectation value for a time constant emission, shown are the 70 % and 90 % probability ranges. 



1.0 




Distance (arcsec) 



Fig. 7. Encircled count fraction perpendicular to the jet axis us- 
ing essentially the photons within stripes 2 and 3 of Fig. |6] The 
MARX simulation shows extended source for comparison. 



woul d be artificially b roader. We adopt a PA of 26 1 ° for the out- 
flow (IPv o et alj|2009l) . 

As we know the absolute position of the jet axis only to about 
0.5", we estimate its position by the mean position of the pho- 
tons perpendicular to the jet axis. The photon centroids along the 
N-S direction diff'er by about 0.14", -0.18" and 0.04" from the 
mean position for the 2001, 2005 and 2009 observation, respec- 
tively, which is well within our estimated astrometric accuracy. 

In order to check whether the X-ray emission region close to 
the driving sources is extended, we show in Figure|7]the distribu- 
tion of the photons perpendicular to the jet axis for the photons 
approximately in stripes 2 and 3 in Fig.|6] Due to the slight incli- 
nation of the jet axis with the x-axis, both regions do not overlap 



exactly. This figure also includes a Marx|f| simulation of an ex- 
tended (1") source for comparison. From this figure it is evident 
that the X-ray emitting region is smaller than 1" and possibly 
smaller than 0.5". We adopt a maximum extent of 0.5" for the 
lateral source extent but note that a source size smaller than this 
value is possible. We also note that the superposition of a number 
of smaller emission regions can mimic the observed photon dis- 
tribution, thus making the physical extent of the emission region 
smaller than 0.5". 

The outer part of the emission region, approximately stripes 
4 and 5, appear extended perpendicular to the jet axis. The lower 
number of counts prevents an estimate of the source size here 
but an extent on the I" level is possible which would imply an 
opening angle of the X-ray emitting jet of about 7°. 

Figure Q shows slight deviations between the individual ex- 
posures, whether these are due statistical fluctuations or due to 
intrinsic changes in the emission region is hard to judge since 
variations are most evident on sub-pixel scales where, even with 
pixel-randomization turned off, the photon locations within the 
individual detector pixel become important. The standard devi- 
ation of the photon distances to the jet axis are 0.65, 0.69, 0.59 
pixel (which is 0.492") for the 2001, 2005 and 2009 observation, 
respectively (energy range 0.5-3.0 keV), i.e., rather similar val- 
ues; the standard deviation for a point source is 0.57. Therefore, 
we will concentrate in the following on changes of the plasma 
properties along the jet axis instead of changes perpendicular to 
it. 

4.4. Knot D in X-rays? 

Between the outer part of the X-ray emission visible in Fig. |6] 
and knot D (or PHK 3), located at a distance of about 13" 
from the driving sources, excess X-ray emission is apparent (see 
Fig. |8]i. However, this excess is not clearly co-aligned with the 
expected postion of knot D. The box shown in Fig. [8] which 

'http : // space . mit . edu/cxc/marx/ 
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Fig. 8. Merged image of the region around HH 154 in the en- 
ergy range 0.4-0.8keV. In order to show the diffuse emission, 
smoothed with a Gaussian kernel of 2.5". The 2005 HST Ho- 
contours are shown. The box connects the inner Ha emission 
with the one of knot D and is of the same size as the box for the 
inner X-ray emission. 



connects knot D with the inner part of the Ha emission, contains 
11 photons in the energy range between 0.5 keV and 3.0keV 
where only 4.5 background photons are expected. With a mean 
energy of 0.92 keV, these photons are also softer than those 
closer to the driving sources. Of these photons, two, three and 
six are from the 2001, 2005 and 2009 observation as expected 
for a constant source due to the different efficiencies. There is 
not discernable spatial evolution. Assuming the absorption of 
knot D also for the region around it {Ay - 2.5 which implies 
Nh - 4.5 X 10^' cm"^), the mean energy indicates a plasma tem- 
perature of about 3 X 10'' K. We caution that merging the obser- 
vations covering about one decade can cancel out any physical 
structure and we regard the above as an indication that X-ray 
emission might be scattered around in this jet on a larger scale. 

Whether the X-ray emission is indeed related to the jet, 
whether it persists over a decade, and what its true space ve- 
locity is cannot be answered from the data. However, we regard 
it unlikely that this X-ray emission represents knots originally 
heated much closer to the driving sources during the last decade 
as space velocities in excess of 1800 km s"' are required to travel 
from the "knot" region (cf., FigjS) to the position of knot D 
within four years. The apparently knot like structure in Fig. [8] 
contains at most two photons from the 2009 observation while 
three are from the previous observations. In any case, the cen- 
ter of this structure is located about 6" from the "knot" region so 
that space velocities in excess of 1000 km s"' are still required to 
travel from the "knot" region to this location. Therefore, this X- 
ray emitting plasms might be the remnant of jet material emitted 
more tha a decade ago. 



5. Discussion 

The interpretation of the X-ray findings depends crucially on the 
required mass loss rate, shock velocity, plasma density and the 
cooling time of the X-ray plasma within the flow. Therefore, we 
start our discussion by deriving estimates for their respective val- 
ues. 



5. 1 . Mass loss rate and shock velocities 

Using the formula given in [Schneider & Schrnit^ (|2008[ eq. 2) 
for the massloss rate required to explain the observed X-ray 
emission we find Mx-ray ~ 1 .5 x 10"'^ Mq yr"' or approximately 
a factor of 10 lower than the inner X-ray emission component 
of DG Tau. This formula assumes that the material observed 
with Chandra is shocked only once, which should be a relatively 
good approximation since high shock velocities are required to 
heat th e material above T > lO^'K. According to Ragaetd] 
( |2002|) . shock velocities of approximately 700 km s ' are needed 
to heat the material to the observed 0.6 keV close to IRS 5, 
while only v shock ~550 kms"' is needed for the lower temper- 
ature at larger distances in case this outer X-ray emitting plasma 
is heated in situ and not the cooled down remnant of the inner 
X-ray component. 

5.2. Densities 

The low densities of the HH flows cannot be directly measured 
with X-rays, however, we can derive a lower limit on the den- 
sity assuming a certain emission volume. From our analysis of 
the source extent in sect. 14.2. H and sect. 14.31 we assume that the 
emitting volume is extended perpendicularly to the jet axis by 
0.5"=:70 AU and along the jet axis by 2" and estimate that at 
least 80 % of the photons originate in this region. These esti- 
mate translates into a maximum volume of 5 x 10"*^ cm"-' for an 
inclination angle of 45°. For a filling factor of unity, the electron 
density of X-ray emitting plasma is then 

/ EM , , 

This value is a lower limit on the density since the plasma might 
be concentrated in individual denser clumps, i.e., the volume fill- 
ing factor could be less than unity. Very close or co-spatial to the 
X-ray emission material of lower temperatures (~ 10"* K) has 
been observed (e.g. [Fe II]), so we consider a filling factor of 
unity rather unlikely since some intermediate temperature mate- 
rial will connect both temperature components. The lower limit 
of the electron density in the outer part of the X-ray emitting jet 
is lower by a factor of about four due to the decreasing surface 
brightness assuming an opening angle of 0° and no change in the 
plasma propertied 

In the optical and near-IR direct density values from line ra- 
tios have been derived mainly for the outer part of the jet where 
densities of a few 10^ cm -' are found, e.g., the [Fe n] Unes im- 
ply Mj, = 7.6 X 10"^ cm"^ dLiseau et a l. 2005). De nsities up to 
10 cm have been derived for the inner 2" bv Itohetal](l2000l) 
from an analysis of [Fe ii] near-infrared lines. Note that the [Fe 
ii] lines have a higher critical density than the [S ii] lines usu- 
ally used for the density measurement. High hydrogen densities 
have also be en found close to the driving sources for a fe w other 
sources (e.g. lMelnikov et al.ll2009t iBacciotti et al.ll2000i) . 

Considering the thermal pressures of these two temperature 
components, we find that the lower limit on the thermal pressure 
of the X-ray emitting plasma {Tx » 7 x 10^ K) is 

Px = 2 nxksTx « 2 x lO"*" dyn/cm^ . (2) 

The high densities of the material observed in [Fe n] results in an 
approximate thermodynamic pressure equilibrium of both com- 

^In sect. 15.41 we argue that the absorbing column density and the 
temperature decrease along the outflow. This changes the lower limit of 
the density in the outer parts of the X-ray emitting jet by a factor of 1.4. 
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ponents, i.e., densities of a few times 10^ cm"-' suffice to pro- 
vide the required pressure at a temperature of T 10'* K. This 
high electron density close to the driving sources supports the 
idea that the density decreases by approximately two orders of 
magnitude within the innermost 5"= 1000 AU (deprojected) as 
estimated from the [Fe ii] and [S ii] lines. 

A conical outflow decreases its density by exactly two or- 
ders of magnitude from 0.5" to 5" for a constant outflow ve- 
locity. However, HH 154 is likely not strictly conical since 
iFridlundet all (120051) noted that the opening angle close to the 
driving source might be as large as 90°, which is consiste nt with 
the estimated lateral jet size of 0.5" bv'P vo etalj (l2002l) at this 
distance. Therefore, the density decrease beyond 0.5" is proba- 
bly less than for a conical outflow as an opening angle of only 
about 3° (see footnoteQ) does not suffice to decrease the density 
sufficiently for large inital opening angles. Thus, it is not clear 
if or where the "cold" jet component is in pressure equilibrium 
with the X-ray emitting plasma. 

Another possibility is that the magnetic pressure supports 
the X-ray emitting volume against expansion. We estimate its 
strength by assuming a plasma-/? of unity {B^ - SnPgas) and 
find B X 6 mG . Such a value is reaso nably close to the driving 
source (see Tab. 1 in Hartigan et al.li2007 ) and requires lower 
densities than pure pressure support. The lower limit would still 
be «e > 4 X 10'*cm"^ if the magnetic field scales with the den- 
sity as expected (B ~ n'' with p - 0.5 ... 1) since measured 
magnetic fields in HH objects indicate 15/iG for n - lOOcm"^ 
dHartigan et al.ll2007h . The interstellar m agnetic field seems too 
weak to colhmate the jet dCabritl |2007|) . however, MHD self- 
collimation is a likely scenario. Therefore, the same wound-up 
(helical) magnetic field which collimates the outflow can lend 
support for the X-ray emitting plasma, largely inhibiting lateral 
expansion of the hot X-ray emitting material. 

5.3. Plasma cooling 

Three processes contribute to the cooling of a plasma: Radiative 
cooling, cooling by expansion and thermal conduction. The pres- 
sure work done by the plasma is SW = pdV and the radiative 
losses are 

6Qrad = -nlVit)A{T)dt, (3) 

where rig is the electron density, V the volume and A(r) is the 
cooling function. The conductive heat flux is given by 

qcond = -K(T) V r , (4) 

where the thermal conductivity according to Spitzer is 

7^5/2 

k(T) = /Co—— erg s"' K"' cm (5) 
In A 

with kq = 1.8 X 10"^ and the Coulomb logarithm In A, which 
describes the collision properties of the plasma and is of order 
10. When the mean free path length for energy exchange is of the 
same order as the thermal scale height, the conduction should be 
approximated by the saturated flux 

qsat = 5(ppcl , (6) 

with Si 0.3 (e.g. iBorkowski et al.|[l989l p is the mass density 
and Cs is the local sound speed). For an estimate of the impor- 
tance of the saturated flux, we assumed a linear temperature de- 
crease. Under these circumstances the Spitzer value exceeds the 

^ The size of the Mach- disk of knot D located at a distance of ~10" 
(0.6". lFridlundetal.l2005h indicates the local size of the jet. 
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Fig. 9. Cooling curves for the different cooling processes. 
Parameters of the models are labeled (density: n, opening angle: 
b, cooling length: [). 



saturated flux on spatial scales of about lOAU for the cooling 
from Ti = 10^ K to To = 10"* K, i.e., the saturated flux should be 
used for these steep gradients (« ^ 10^). 

Thermodynamics states that the energy change of a plasma 
cell is described by 

dU + 6W^6Q^ SQrad + qcond ■ A dt (!) 

with the internal energy U - aNkT (a = 3/2 for a fully ionized 
plasma), the particle number A^, the Boltzmann constant k, the 
temperature T and A is the surface area through which heat con- 
duction proceeds. We us e p - IngkT in the expression for the 
pressure work and follow iGiidel et al.l (l2008l) by writing eq.|7]as 

a + — --{ H -A WT\dt, (8) 

T(t) V \2kT(t) IrieVkT In A / 

where we used A^^. = «p • V and note that this expression holds 
only in the presence of sufficiently small temperature gradients. 

In order to estimate the relative importance of the three cool- 
ing terms, additional information is needed, in particular, the 
opening angle of the X-ray emitting jet, its density structure, the 
temperature gradient, the surface for the heat conduction, which 
would include the magnetic topology and the properties of the 
environment, e.g., its ionization. These quantities are not avail- 
able for the X-ray emitting part of the jet. We therefore decided 
to give some order of magnitude estimates for the cooling times 
of the individual processes ignoring contributions of the other 
ones. As we will see, there are distinct regions in the parameter 
space where each process seems to dominate, so we regard this 
approach reasonable. 

Figure |9] shows the cooling curves for the different processes 
assuming different parameters for the jet. For radiative and con- 
ductive cooling, the mapping of time to distance in this figure 
depends on the actual, deprojected space velocity of the plasma. 
A rough estimate is 0.3" yr"', which implies that today's inner 
X-ray emission will reach the position of today's outer emission 
in 15 years. Adiabatic cooling, on the other hand, does not de- 
pend of the outflow velocity but only on the initial cross-section 
of the plasma and on the opening angle. 
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5.3.1. Adiabatic cooling 

Protostellar jets usually show an approximately conical struc- 
ture at some distance from the driving sources so that the flow 
expands mainly perpendicular to the jet axis. In the limit of adi- 
abatic cooling 

r V"'"i = const (9) 

holds. Since we do not observe local temperatures, we have to 
average the temperature, weighted by densitjj^, over the volume 
used to measure the temperature. We use the following approxi- 
mation for the volume of the plasma cell 

V ^nlin^ + nanbf , (10) 

where r(f) = v ■ f is the position along the jet axis measured 
from the initial distance, while ro = 0.25" is the initial jet radius 
at this position and 2 ■ b the opening angle (I is the length of 
the cell along the jet axis). The initial cross-section is fixed and 
the temperature decrease depends only on the position along the 
outflow. From the size of the Mach disk about 10" from the 
driving sources, we estimate an opening angle of 3 - 10° for the 
flow, where the separation of the working surface and the Mach 
disk argues for values closer to 3°. Diff'erent outflow velocities 
would change the curve for the expansion cooling in Fig. |9]but 
would not lead to another spatial temperature structure, because 
the dependence on v cancels out in the eq uations. 

As described by Giide l et al.l (l2008l) . the expansion addi- 
tionally reduces the density of the emitting plasma and thereby 
lowers the number of emitted photons more strongly than ex- 
pected on the basis of the temperature decrease alone. For a con- 
sistency check, we calculated the expected number of photons 
at 3.5" from the driving sources from the ratio of the emis- 
sion measures at 0.5" and 3.5" and the drop in temperature. 
Assuming constant absorption, we expect a drop in photon num- 
ber by approximately a factor of about 6 from 0.5" to 3.5" for 
an opening angle of 3°, which is approximately compatible with 
the observed value. The larger opening angle of 10° would re- 
duce the photon number more strongly, i.e., the combination of 
the temperature and density decrease reduces the expected pho- 
ton number by about 200 for the same distance. 

5.3.2. Radiation cooling 

We solved eq.lHusing th e cooling function of Chianti version 6.0 
dPere et al I1997II2OO9I) assuming half solar metallicity. Figure|9] 
shows two cooling curves for radiative cooling. According to 
eq. [8] the cooling time depends linearly on the density. It is 
clear that radiative cooling does not contribute significantly to 
the cooling as long as the density does not exceed n ^ 10* cm"^. 

5.3.3. Conductive cooling 

Magnetic fields are essential for the launching of jets, but even 
at greater distances, small magnetic fields (~ lOOyuG) influence 
the jet dynamics (Hartiganet al. 2007). They can also strongly 
suppress heat conduction perpendicular t o the field lin es even 
for weak fields (~ lju G, see eq. 5-53 in ISpitzeil[l962l) . In the 
presence of turbulent magnetic fields, heat conduction might be 
suppressed by about two orders of magnitude or even enhanced 
relative to the Spitzer value (e.g. Naravan & Medvedev 2001; 
ICho et al.l 120031 : lLazarianll2006l) depending on the scale of the 

^Note that EM = rP-V = nN with a constant number of particles N 
in each cell. 



turbulence. We regard it as plausible that heat conduction works 
most efficiently along the jet axis while it is suppressed by some 
kind of magnetic field perpendicular to the jet axis. The Spitzer 
value for the heat conduction assumes an ionized plasma, which 
might not be entirely true throughout the jet, however, a consid- 
erable amount of ionized material should be present close to the 
X-ray emitting plasma. Given these uncertainties, we estimate 
conductive cooling by 




given in Orla ndo et al ] 12005). We show in Fig.|9]a cooling curve 
by numerically integrating the conductive cooling for a fixed 
density n and for cylindrical geometry {V ~ A ■ I). The effect 
of the conductive cooling depends on the density of the plasma 
and on the temperature gradient, i.e., on the cooling length (we 
assumed 600 AU for a temperature decrease from 0.6 keV to 
0. 1 keV). The curve shown in Fig. |9]is intended to give a rough 
impression of this effect and we caution that the provided esti- 
mate for the conductive cooling might be off by orders of mag- 
nitude in some scenarios, e.g., for turbulent magnetic fields. 

5.3.4. Cool conclusions 

From Fig. |9]it is clear that cooling by expansion dominates over 
radiation. Whether conduction is important depends on the den- 
sity, the temperature gradient and the magnetic field configura- 
tion. When no heat is transferred perpendicular to thejet axis, we 
expect adiabatic cooling to dominate. We will therefore focus on 
that cooling process in the following. 

5.4. Trend in mean energy 

Judging from the absorption value (Nh ~ 1.4 x 10^^ cm"^) close 
to the driving sources, th e visual ab sorbing magnitude is approx- 
imately Al/ 8 mag dVuong et al.ll2003i) . Iltoh et alj (1200 0) an- 
alyzed the [Fe 11] 1.644 jum/1.257 jum line-ratio as an estimate 
for the evolution of the visual extinction along the jet axis and 
found a value of Ay » 7 for distances greater than 1 " from IRS 5. 
This is compatible with our estimate from the X-ray spectrum. 
Therefore, we expect that the coiTelation of Ay and Nh holds 
in the jet region. Closer to the driving sources the extinction in- 
creases up to Ay «i 2 1 . iFridlund et al.l (l2005h also d erived Ay- 
values, which are slightly lower than the values of Ito h et al.l 
(2000) for the optical knots located further downstream. They 
find that the extinction decreases slowly towards the outer knots 
where Ay « 2 - 3 is found. At the position of knot F which 
is closest to the driving sources, they estimate an absorption of 

Ay > 4. 

Assuming that the absorbing column density decreases from 
Nh = 1.4 X 10^^ cm"^ at 1" from the driving sources to Nh - 
8 X 10^' cm"^ at 5", the plasma temperature still has to decrease 
from 0.7 keV to 0.4 keV in order to explain the decrease in the 
mean energy. We therefore conclude that the temperature de- 
crease of the X-ray plasma significantly contributes to the soft- 
ening of the photons beyond a: 1". It is possible that some emis- 
sion is coming from the innermost 1 " where it is more strongly 
absorbed. For an upper limit on the temperature change along 
the flow axis, we fix the absorption to Nh = 8 x 10^' cm"^, 
which requires, according to Fig. |5] that the plasma tempera- 
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ture decreases from approximately 1.7 keV to below 0.2 keV 
within 3". 

An opening angle of a few degrees reduces the plasma tem- 
perature along the outflow as required by the above estimates, 
consistent with the assumption that adiabatic cooling dominates 
the plasma cooling. As discussed in sect. 15.3. li the reduction of 
photons along the jet axis is compatible with a shallow opening 
angle. 

As we have no good estimate for the cooling time of the 
plasma, it is possible that cooling is very efficient and re-heating 
of the X-ray emitting plasma along the inner few arcsec is re- 
quired in order to produce observed X-ray emission at larger 
distances from the driving sources. However, the temperature 
decrease along the jet axis remains virtually constant over 10 
years of observation so that also the hypothetical re-heating must 
be relatively constant over this period of time. Large individual 
blobs with largely varying velocities, ejected every few years, 
would probably produce a more variable temperature structure. 
Smaller, unresolved internal shocks could, on the other hand, be 
present so that the cooling time of the plasma can be shorter than 
derived from the decrease in the mean energy. 

5.5. The inner emission component 

The most striking feature observed in all observations is the 
existence of a luminous X-ray emission region close to the 
driving source(s). The peak of this feature is approximately 
0.5-1.0" offset fromL 1551 IRS 5 (R.A.(2000)=04h3 1^34^5, 
Decl.=18°08'05'.'04) towards the south-west. Its luminosity and 
temperature remains virtually constant within a timespan of 
about nine years. 

Therefore, it seems necessary to review the arguments which 
prohibit an association of this inner component with one or both 
of the central driving sources themselves. Essentially, these are 
(a) the astrometry and (b) the absorption. Concerning (a), the 
centroid of the inner component is placed at least 0.5" from the 
location of radio sources in every of the three available Chandra 
observations. Although this is marginally compatible with our 
estimated astrometric accuracy, this is unlikely to be caused by 
a repeated incorrect pointing of the satellite since the centroids 
(using the photons in the inner r = 1" circle) of the inner com- 
ponent match to within 0.3" for all observations. As to (b), the 
interpretation of the scattered light and the non-detection of di- 
rect emission places a firm lower limit of Ay » 20 on the ab- 
sorption towards IRS 5 (IStocke et al.lll988T) . This translates into 
an absorbing column density of Nh » 4 x 10^^ cm"^, at least 
three times higher than measured. The measured value, on the 
other hand, corresponds well to the one estimated for the inner 
jet from the near-infrared line-ratios and suggests that X-ray and 
ne ar-infrared emission spatially coincide. 

iBallv et all (l2003h sketched a possible scenario in which 
scattered X-rays are responsible for the observed X-ray emis- 
sion. However, these authors concluded that this option is 
less likely and we agree with their evaluation. The parame- 
ters required for this scattering scenario, e.g., densities of « ~ 
10^ cm"-*, are not strictly ruled out, but would be exceptional for 
protostellar jets. Furthermore, the trend in the mean energy also 
argues against the scattering scenario, because Thomson scatter- 
ing is not very sensitive to the scattering angle and independent 
of the wavelength. Dust scattering, on the other hand, is by far 
not sufficient to explain the observations using the usual conver- 
sion factors. 

We therefore associate this feature with th e apparen t ly sta- 
tionary [Fe ii] emission complex observed by iPvo et al] (l2009l 



l2005Ll2002h . Their [Fe n]Al .644 //m data, obtained over a times- 
pan of four years, show an apparently stationary component 
close to the driving sources. Pvo et al. (2009) already proposed 
that the inner X-ray emission component is associated with the 
innermost [Fe ii] emission peak, called PHKl (distance to IRS 5: 
1.1"). The total flux in [Fe ii] (~ 5 x 10"'^ erg/s/cm^) is within 
an order of magnitude comparable with the X-ray flux (0.5- 
lO.OkeV: 4x10"'^ erg/s/cm^ orunabsorbed3xlO"''* erg/s/cm^). 

The [Fe ii] emission at PHKl is dominated by low veloc- 
ity material with v ~ -60- ■ ■ - 150 km/s (or deprojected v = 
-85 . . .212 km/s, for an inclination of / a 45°). Interestingly, 
the post-shock velocity of a shock with an initial velocity of 
v=700 km/s is 17 5 km/s, i.e., within th e range of the low ve- 
locity component. lPvo et al.l (I2002il2009l) noted that the velocity 
dispersion of the low velocity [Fe ii] emission decreases with in- 
creasing distance to the driving source, which they interpret as a 
collimation of the outflow. The large opening angle of the flow 
close to the driving sources and the shallow opening angle fur- 
ther downstream support this interpretation. Collimation might 
also be responsible for the X-ray production and would natu- 
rally explain their stationary appearance. Since the high- velocity 
[Fe ii] material appears approximately where the X-ray emis- 
sion disappears, it is tempting to associate this material with 
outflowing plasma not as strongly shocked as the X-ray emit- 
ting material. However, the total mass-loss derived for the X-ray 
component is lower than for the optical part of the jet so that 
it remains unclear whether the "absence" of highly blueshifted 
emission close to the driving sources is somehow connected to 
the existence of X-ray material, i.e., if a large fraction of shocked 
high-velocity material reaches X-ray emitting temperatures. 

5.6. The extended or outer emission component 

During the 2005 ACIS-I observation, an enhancement of pho- 
tons two arcseconds downstream from the bulk of the X- 
ray emission is evident (Favata et al.l I2006h . We estimated in 
sect. 14.2.21 that the low number of counts in the corresponding 
region might be a statistical fluctuation. Nevertheless, it is still 
possible, if not even physically plausible, that the elongation dif- 
ferences are caused by a transient X-ray emitting knot, possibly 
comparable with other X-ray emitting knots within HH objects. 

Concerning the position of this blob, its 2005 position co- 
incides with one or all emission peaks in the "F"-complex 
(iBonito et al.l 12008). The space velocities measured in this re- 
gion range from ~ 100 kms~' to 500 k ms~^ in optical forbid- 
den emission lines dFridlund et al.ll2005l) . At this distance from 
the driving source the high-velocity component in [Fe ii] be- 
comes dominant over the low- velocity compo nent, and veloc - 
ities up to 5 00 kms ' have been measur ed (iPvo et al.ll2005h . 
Furthermore, iPvo et al.l (l2005l) noted that the outer high veloc- 
ity [Fe ii] component might exhibit time-variability at approx- 
imately the same time of the appearance of the X-ray knot. 
However, the F-complex did not change mu ch during this time 
in its optical appearence jBonito et al.ll2008]) . 

The distance traveled by this hypothetical knot between 
2005 and 2009 would be approximately 0.6-3.0" (100 kms"' 
. . . 500 km s"'). In the corresponding regions zero or one photon 
are recorded during the 2009 ACIS-S exposure, which is more 
sensitive than the previous ACIS-I observation at low energies. A 
maximum of seven photons can be attributed to the 2005 "knot" 
(cf. Fig. |2] where the most favorable geometry is sketched), so 
that the luminosity of this knot must have decreased during the 
last four years if it indeed existed. For the interpretation of these 
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phenomena, the cooling time of the X-ray emitting plasma is 
crucial but unfortunately not known with the required precision. 

When the cooling time of the hot plasma is short compared 
to the travel time to the outer locations of X-ray emission, it is 
impossible that the material is only heated close to the driving 
sources and then just cools while it is flowing outwards. Internal 
shocks are a natural explanation for the re-heating, but the ob- 
servations require a nearly constant decrease in shock velocity 
with increasing distance to the driving sources (see the trend in 
the mean energy, sect. 15.4b . 

If, on the other hand, the plasma is not significantly re-heated 
while flowing outwards (no internal shocks), either a variable 
mass outflow or a statistical fluctuation are responsible for the 
apparent knot. In any case, a rather constant temperature close 
to the driving sources is required, which translates to a constant 
shock velocity for shock heating. In both cases, a shallow open- 
ing angle is mandatory and magnetic fields probably suppress 
heat conduction efficiently. The decrease in plasma temperature 
along the flow reflects the cooling time of the plasma and ex- 
plains, why no emission is detected at larger distances from the 
driving sources. 

5. 7. Comparison with DG Tau and otiner jet X-ray sources 

The so-called TAX sources (Giidel et al. 2007b), which show 
spectra composed of two emission components with vastly dif- 
ferent absorbing column densities, e.g. from an embedded star 
and a less embedded jet, show a striking similarity to HH 154. 
We associate the beehive proplyd and similar COUP sources 
with this group (K astner et alJl2005h . 

The comparison of HH 154 with the X-ra y emitting jet of 
DG T au (Gudel et al. 2005, 2007b, 2008; Sch neider & SchmittI 
120081: iGUnther et al.ll2009i: iGudel et al.ii201 1) is particularly in- 
teresting as multi-epoch high spatial resolution X-ray observa- 
tions are also available. The important similarity between both 
jets is that the majority of the X-ray emitting plasma is located 
close to the driving source(s) in all observations. We will there- 
fore focus our comparison on this innermost part of the outflow, 
the resolved outer part of DG Tau's jet is probably more re- 
lated to the outer part of HH 154 (knot D, see sect. 14.4b . The 
inner part of the outflow also appears qualitatively similar i n 
the [Fe II] position-velocity maps (e.g. iPvo et"ani2003L l2009l) . 
Both outflows exhibit a low velocity component close to the 
star, while the high velocity component is located further down- 
stream. Other jet-driving sources like HL Tauri and RW Au rigae 
show a slightly different pattern in [Fe II] (iPvo et al.ll2006 ) with 
the low velocity component at a larger distance from the stellar 
position. 

The cooling times of the X-ray emitting plasma, however, 
appear to differ between the two jets. While we estimated in 
sect. l5.4l that the temperature gradient along the outflow indicates 
a cooling distance of the X-ray emitting plasma in the range of 
a few lO OAU., the essentially unresolve d inner component of 
DG Tau dSchneider & Schmittll2'008; .Gudel et al.ll201 1} implies 
a much shorter cooling time for DG Tau. This might be related 
to the higher minimum electron density of the X-ray emitting 
plasma of this component in DG Tau (n^ > 10^ cm"^) compared 
to HH 154 (rie > 1.2 x 10^ cm"^). While the X-ray emitting 
plasma of HH 154 can be in pressure equilibrium with the cooler 
jet material, this seems much less likely for the X-ray emitting 
jet of DG Tau. The shorter distance of the inner X-ray jet com- 
ponent of DG Tau to its driving source compared to HH 154 
might be the cause for the higher density of DG Tau's jet X- 
ray emission. Assuming that the heating happens more closely 



to the driving source for DG Tau's jet, a large opening angle of 
the outflow can explain the essentially point-like appearence of 
the X-ray emission in DG Tau as cooling by expansion would 
reduce the plasma temperature sufficiently fast. The variable X- 
ray luminosity without measur able differences in the position of 
the inner X-ray jet component (iGiidel et al.l201 ll) also argues for 
a short cooling distance of this X-ray component as otherwise at 
least shifts in the mean position of the X-ray emission should be 
measured assuming the proper motion of the optically observed 
knots (~ 0.3" yr For HH 154 a shallower opening results in 
a larger extension of the X-ray emission. One can speculate that 
the X-ray emitting plasma contributes to the expansion of the jet 
close to the driving source in DG Tau while that is not as clear 
forHH 154. 

An obvious difference between the X-ray jets of DG Tau and 
HH 154 is the different plasma temperature close to the driv- 
ing source. However, the soft component of the beehive pro- 
plyd has approximat ely the same plasma temperature as HH 154 
(lKastneretal.ll2005h . For DG Tau, the inner X-ray component is 
only marginally hotter than the outer, resolved one at a distance 
of 5" from DG Tau, while for HH 154 the outer parts of the 
jet (e.g. knot D or the outer stripes in Fig. |6) seem to be cooler 
by a factor of two compared to the emission close to the driv- 
ing source. Since adiabatic cooling decreases the temperature 
strongly along the flow, it seems unlikely that the same plasma 
found close to the driving source is also responsible for the outer 
X-ray component of the DG Tau's jet or knot D in HH 154, and 
additional internal shocks cause the high temperatur es further 
out wards (e.g. as in the models of lBonito et al.ll2010bh . 

IStelzer et al.l (l2009h recently detected the appearance of an 
X-ray knot after an FU Ori like outburst of Z CMa, i.e.. X-ray 
emission located about 2000 AU from the driving source. This 
emission is much farther out than in HH 154, which can either 
indicate that the lifetime of such knots might be relatively long 
or that strong shocks also occur further outwards in the flow. 
As in Z CMa the X-ray emission in other HH objects is located 
further out along the jet. If the correlation of the X-ray emission 
region with one of the working surfaces in these outflows is an 
observational bias, such that X-rays between the knots are less 
likely to be recognized as associated with the HH object or if 
it is an intrinsic feature of the X-ray production mechanism is 
currently not clear. 

In summary, two properties make X-rays from the proto- 
stellar jets identifiable, (a) a TAX-like spectrum and (b) X-ray 
emission displaced with respect to the driving source. Strictly 
speaking, HH 154 belongs to the TAX class of objects and the 
strong absorption of the driving sources makes the detection of 
extended emission possible. We speculate that the X-ray mor- 
phology of HH 154 also applies to the other TAX sources where 
the angular resolution is insufficient to resolve all details, e.g. 
due "contamination" by stellar emission. As the outflow rate re- 
quired to produce the observed X-rays in these objects is lower 
than estimated from the optical emission by a few orders of mag- 
nitude, it seems likely that such a high velocity component has 
escaped detection in other wavelength regimes, but is still pow- 
erful enough to lead to the observed X-ray fluxes. 

6. Model implications 

After the discover y of X-ray emission offs et from the driving 
sources (IRS 5) bv lFavata et aP (12002 ') and'Ba llv et aP (l2003h . 
a variety of models have b een proposed in or der to explain this 
phenomenon. In particular iBallv et al] (l2003b described an en- 
semble of possible explanations covering a broad range of pos- 
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sibilities including, e.g.. X-rays from the driving source (s) re- 
flected into the line of sight from the outflow cavitv. ^onito et all 
(2010a) performed detailed magneto-hydrodynamical simula- 
tions of a jet with a variable outflow velocity focusing on high 
velocity shocks within the outflow and the associated X-ray 
emission. These authors discuss four different scenarios for the 
evolution of the X-ray emission of HH 154, one of which con- 
siders a stationary source. 

We concentrate on models related to the apparently station- 
ary X-ray emission complex since the majority of the X-ray 
emission is related to this complex. The scenarios related to the 
interpretation of the outer or extended X-ray component were 
discussed in sect. l5.6l and depend crucially on the unknown cool- 
ing time of the X-ray emitting plasma. 

6.1. A jet with random ejection velocity/Internal shocks 

This is the model discussed by iBonito et alJ (l2010allbl) . Their 
simulations can produce an emission complex close to the driv- 
ing source, when a recently ejected faster blob overtakes a 
more slowly moving blob (cf. their Fig. 2). In the absence of 
strong cooling, the proper-motion of such a knot would be de- 
tectable with the available high resolution X-ray observations. 
Concerning the inner, apparently stationary source, the models 
of|Bonito et al. (2010b) predict that the most probable position 
of a shock is close to the driving sources. 

The virtually constant X-ray luminosity and the relation to 
the constant [Fe II] emission argues against a strong variation 
of the shock velocity or location. One solution for these dis- 
crepancies is a relatively regularly modulated jet so that a con- 
stant luminosity might be mimicked by the superposition of a 
roughly constant number of smaller shocks formed close to the 
driving source. The trend in the mean energy would then reflect 
the cooling of these smaller individual blobs, while they travel 
alo ng the outflow (e. g., see also the sub-radial blobs modeled 
by lYirak et ai]|2009l) . Another possibility is that variations in 
the shock properties are hidden by the low photon numbers and 
the inner X-ray emission is caused by larger knots shocked close 
to the driving source. The constant appearance would then not 
reflect a constant outflow but would rather be a chance coinci- 
dence. 

In this model the absence of X-ray emission farther down- 
stream would be explained either by a lower density or a low 
temperature of the plasma inhibiting its detection. The opening 
angle of the optical jet is roughly consistent with this picture. It 
requires, however, that strong shocks at larger distances from the 
driving source are less probable t han close to t he driv ing sources, 
which is true for the models of Bonito et al. In case of 

a short cooling time the decrease of the plasma temperature re- 
flects a decreasing shock velocity with increasing distance to the 
driving sources, which would provide another explanation for 
the non-detection of X-ray emission farther downstream. In any 
case, the observations clearly show that the heating to X-ray tem- 
peratures is a function of the position along the flow. 

6.2. Base/collimation shock 

Guided by the first Chandra observation. lBallv et al.l (l2003h pro- 
posed that some kind of stationary base-shock can explain the 
observed X-rays, either independently for each driving source 
or at the envelope of both sources. In these scenarios, the mag- 
netic fields can coUimate the outflow and can also support the jet 



against the thermodynamic pressure of the hot X-ray emitting 
plasma. 

This scenario requires lower velocities than the internal 
shock model, but still higher than detected in available spectra. 
The "deflection" angle might be relatively large (~ 45°) as the 
opening angle close to the driving sources might also be large so 
that flow velocities of 10^ kms"' suffice for the X-ray produc- 
tion. Also, the concentration of the X-rays within a rather small 
volume close to the driving sources and the virtually constant 
X-ray luminosity are a natural consequence of this scenario . The 
base-shock scenario does not inhibit jet mass flux variations and 
is consistent with the observations as long as the amplitude of 
these variations is small enough. For sufficient dumpiness and 
blob ejection cadence, the base-shock model and the internal 
shock model become indistinguishable and share the possibility 
for small amplitude time variability. 

Concerning the location of the base-shock, it seems pleas- 
ing to attribute the brightest X-ray spot to the locat ion of the 
base-shock consistent with the [Fe ii] observations of iPvo et"an 
(2009). However, the increasing absorbing column, which seems 
to cause the hardening of the photons in the innermost part of the 
flow, can absorb the soft X-rays closer to the driving sources. 
Thus, it is possible that the true location of the shock region is 
hidden behind a larger absorbing column and located closer to 
the driving source. 

The trend in the mean energy is also a natural outcome of this 
scenario, when the plasma is heated to X-ray emitting tempera- 
tures close to the driving sources and cools while flowing out- 
wards. Adiabatic cooling, providing an upper limit on the cool- 
ing time, is approximately consistent with the observed trend in 
the mean energy. 

6.3. Processing jet 

Jet precession seems to be required for some of the observed 
jets (e.g. HH 34 lMasciadri et alJ2 002). The precession times are 
usually rather long (> 10^ yrs), therefore, the change in outflow 
direction for HH 154 would be small between the 2001 and the 
2009 Chandra observations. Still, some kind of a drilling ef- 
fect might be present. A constant flow hitting different parts of 
the envelope would lead to a constant appearance of the inner 
emission component. As the opening angle (~ 90°) close to the 
source is probably large compared to the expected precession 
angle, we regard it as less likely that the direction change of the 
outflow is responsible for the X-ray emission. 

6.4. Stellar wind 

The solar wind has roughly the temperature and velocity ob- 
served for HH 154's X-ray emission. We can imagine that 
during the early stellar evolution the outflow rates of stellar 
winds are much higher than for the present-day Sun and that 
the same process leading to the coUimation of the slower out- 
flow components collimates the stellar wind. A stellar wind 
might be important for the angular momentum problem, but 
cannot be responsible for outflow rates above 10'*^ Mp/ year 
due to the resulting e xcessive X-ray emission ( Decamphlil98Tl 
iMatt & Pudritd l2b08h . However, the outflow-rate required for 
the observed X-ray emission i s orde rs of magnitudes lower and 
iGomez de Castro & Verdugol (l2007l) found evidence for a stel- 
lar driven wind for RY Tau. These authors suggest that the su- 
perposition of many individual small-scale outflows from the 
stellar surface leads to observed morphology of the FUV lines. 
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Therefore, a stellar wind, while not responsible for the main out- 
flow, might provide the required high temperature plasma close 
to the driving source. This stellar wind would not require high 
outflow velocities for the shock heating since it is already of ap- 
proximately the correct temperature when launched. However, 
tile association of X-ray emission with shocked material at other 
wavelengths makes this explanation less likely but would be an- 
other possibility explaining the constant appearance. 

7. Summary 

Our new, third epoch Chandra observation clearly shows that 
the process responsible for the X-ray emission in HH 154 is con- 
stant over at least one decade. The position, the luminosity and 
temperature of the X-ray emission are virtually the same in all 
observations. Whether differences between the observations are 
statistical fluctuations or intrinsic differences in the flow cannot 
be definitely decided due to the low count statistics. From the 
trend of the mean energy along the jet axis, we show that the 
plasma is cooler at larger distances from the driving source. 

We discus several models and find that a standing shock 
most naturally explains tiie observed morphology given the con- 
stant total X-ray luminosity. The location of the X-ray emis- 
sion, where the outflow is likely collimated, and its stationary 
appearance argue for this model. Depending on the details of the 
plasma cooling, the trend in the mean energy can be naturally 
explained in this model. The features of the X-ray emission can 
also be explained in terms of a pulsed jet, where internal shocks 
cause an apparent stationary X-ray source as the most proba- 
ble location of an X-ray emitting shock is close to the driving 
source. The trend in the mean energy might then reflect lower 
shock velocities or the cooling of tiie plasma depending on the 
detailed cooling times of the X-ray emitting plasma. The ex- 
istence of knots within protostellar jets is usually attributed to 
time-variable outflows, therefore, such a model is attractive, but 
it requires a rather regularly modulated flow, since the position, 
the temperature and the luminosity appear constant. Variability 
at larger distances from the driving sources might be present and 
can be explained either by local shocks or variations of the mass 
loss rate. 

A comparison of our new results for HH 154 with other X- 
ray emitting jets, in particular with DG Tau, the only nearby 
jet X-ray source where multi-epoch observations are available, 
shows that soft X-ray photons close to the driving source are not 
unique to HH 154. Therefore, the necessary heating apparently 
takes place very close to the driving source within the outflow. 
With an increasing number of X-ray observations it becomes in- 
creasingly clear that the origin of the X-rays is tightly connected 
to the flow properties within the innermost few 10 AU, either due 
to inhomogeneities in the outflow or by the coUimation process. 

New sensitive X-ray observations of HH 154 with a higher 
cadence are required to decide whether variations on shorter time 
scales are present and could therefore discriminate between tiie 
base- shock and the internal shock model. 
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